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Using a combination of the variation approximation (VA) and direct simulations, we consider the 
hght transmission in nonlinearly amplified bulk media, taking into account the localization of the 
gain, i.e., the linear loss shaped as a parabolic function of the transverse radius, with a minimum at 
the center. The balance of the transverse diffraction, self-focusing, gain, and the inhomogeneous loss 
provide for the hitherto elusive stabilization of vortex solitons in a large zone of the parameter space. 
Adjacent to it, stability domains are found for several novel kinds of localized vortices, including 
spinning elliptically shaped ones, eccentric elliptic vortices which feature double rotation, spinning 
crescents, and breathing vortices. 
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Introduction and the model. Self-collimated and self- 
guided light beams organized as spatial solitons in bulk 
media are subjects of great interest [l[. In addition to 
their significance for fundamental studies, spatial solitons 
can find applications to the design of all-optical data- 
processing schemes. In particular, laser systems may run 
on dissipative optical solitons, which are described by 
the complex Ginzburg-Landau equations (CGLEs), typ- 
ically with the cubic-quintic (CQ) nonlinearity, taking 
into account the saturable nonlinear gain [2]. Crucially 
important to the applications is the stability of dissipa- 
tive solitons. Families of stable spatial solitons with two 
transverse dimensions (2D) exist in materials character- 
ized by a saturable nonlinearity that compensates the 
diffraction [l|, y] . The additional condition of the balance 
between losses and gain reduces the family to isolated so- 
lutions, one of which may be stable as an attractor [J, |5|. 

The complexity of the CGLEs does not admit exact 
solutions, with rare exceptions [a, Q- Nevertheless, an 
analytical approximation, which provides clues to find- 
ing dissipative solitons in the numerical form, has been 
developed in the form of the variational approach (VA) 
adopted for dissipative systems Q- Parameter domains 
where 2D and 3D solitons are stable in the CQ-CGLE 
and related models have been outlined by means of these 
methods 0,il. 

Objects of fundamental significance which remained 
elusive in the studies of dissipative solitons are solitary 
vortices. In previous works, such vortices were found 
to be stable only in the presence of the diffusion, which 
occurs in other physical media, but does not appear in 
optical systems [^, [Tl[ , nor in dissipative models of Bose- 
Einstein condensates (BECs) 0, [l3|- The stabilization 
of multi-peak patterns carrying vortical phase patterns 
was also demonstrated by periodic potentials that may 
be created in the laser cavity @- Nevertheless, the most 
fundamental circular ( "crater-shaped" ) localized vortices 



were not found yet as stable objects in models of optical 
media. In this work, we report (for the first time, to our 
knowledge) stable fundamental vortical modes in a phys- 
ically relevant model where the stabilization is provided 
by an axisymmetric modulation of the linear loss, with a 
minimum at the center. In the experiment, the loss mod- 
ulation may be induced by a localized gain which partly 
compensates the uniform loss (an "iceberg of gain" sub- 
merged into the "sea of loss"). In fact, the gain applied to 
a laser cavity is always localized (in particular, due to the 
finite width of the pumping beam). Besides the optical 
media, similar modulated loss profiles can be engineered 
in BEG, where they help to support various matter-wave 
modes Q. 

We also demonstrate that, when the axisymmetric vor- 
tices lose their stability, they give rise to other novel types 
of vorticity-carrying modes. These are spinning ellipti- 
cally deformed vortices, eccentric elliptic vortices which 
feature spinning and precessing, revolving crescents, and 
breathing vortices. Thus, the use of the 2D "iceberg- 
shaped" gain profile opens the way to the previously in- 
accessible class of stable vortical modes in bulk optical 
media. 

We adopt the (2+l)D form of the CQ-CGLE govern- 
ing the evolution of the optical field in the bulk medium, 
E (z, X, y), along axis z, and the diffraction in the trans- 
verse plane @, 0] : 

lE, + i {E,, + Eyy) + (1 - le) \E\^E ~ {v ~ ly.) \EtE 

(1) 



= —ig{r)E, r = y^x^ + y^, 



where positive coefficients v, e and fi account for the sat- 
uration of the Kerr nonlinearity, cubic gain, and quintic 
loss, respectively. As said above, the linear-loss coeffi- 
cient is modulated along the radial coordinate, g(r) — 
7 -I- Vr"^, with 7, y > 0. Note that the limit form of a 



similarly organized ID model, with the linear gain con- 
centrated at a single "hot spot" (approximated by the 
delta-function) and cubic nonlinearity, admits exact sta- 
ble solutions for pinned solitons [6]. 

The variation approximation. First, we aim to obtain 
approximate analytical results, using the VA technique 
[3|. To this end, we adopt the ansatz for vortex soliton 
with topological charge 1, 
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where 9 is the angular coordinate, normalizing factors 
are A^ = 3/(2-\/2) and i?* — 16/9, while variation pa- 
rameters are amplitude A, radius R, radial chirp C, and 
phase ^. The total power corresponding to ansatz ^ 
is P = 27r/g°° \E{r)\'^rdr = P.A^R^, with P, = ttAIrI 
Skipping straightforward details, the VA leads to the sys- 
tem of evolution equations: 



dA A /5 



dz i?? V2 



eA" - 2/xA^ -2C]-A{-f + RiVR'), (3) 



f = |(^^-|-^'-^'-«2™^ 



dC _ 1 / 1 A^ 

H^ ~ ~m \2m^ 2m 



2m 



2C' 



(4) 



(5) 



Steady-state solutions are obtained as fixed points (FPs) 
of Eqs. ©-(H]). In the case of small chirp, it follows from 
Eq. ([5|) that the width of the stationary state depends 
only on its amplitude, R^ = {A^ — vA^)'^, as in con- 
servative systems. In the first approximation, the small 
chirp, which makes the dissipative solitons different from 
their conservative counterparts, is obtained from Eq. Q: 
C = A^{e - piA^)l^ -K RiVl [2^2 (l - vA^)\ . 

Solving the FP equations following from Eqs. ([3])- 
([5]) analytically, we find two different solutions for A, as 
shown in the inset to Fig. 1. As per general principles 
of the analysis of dissipative systems, the solution with 
larger A may be stable, while its counterpart with smaller 
A corresponds to a separatrix between the two attractors 
Q - the stable vortex and trivial state, E = Q. The solu- 
tion with larger A is characterized by C < 0, which also 
is a necessary stability condition for the steady state [3] . 
The stability of the FPs was checked through the com- 
putation of eigenvalues for small perturbations, see the 
results in Fig. 1. 

Numerical results. The full stability area for the vortex 
solitons was identified by means of direct simulations of 
Eq. ([T]), using the VA-produced ansatz as the input. The 
VA predicts the qualitative shape of the area correctly, 
although overestimating its size. A typical example of 
the fast formation of a stable vortex (which completes 
by z = 20, i.e., after passing < 4 diffraction lengths) is 
displayed in Fig. 2. 
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FIG. 1: (Color online) Stability regions of vortex solitons in 
plane (e, V) of the loss-modulation parameters. In this figure 
and below, other coefficients in Eq. ((T| are 7 = 0.29, /x = 1.4, 
V = 0.4, which adequately represent the generic situation. 
The VA predicts the stability in the shaded (yellow) area, 
and the corresponding amplitude displayed in the inset. 
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FIG. 2: (Color online) The self-trapping of a stable vortex 
from the input predicted by the VA, at e = 2.325 and V = 
0.03. 



In the "filamentation" region in Fig. 1, the modula- 
tional instability breaks the vortex into two segments, at 
z ~ 1000 (closer to the stability border, this distance in- 
creases to z ~ 3000). This outcome of the evolution is 
explained by the fact that, in this region, the total power 
of the vortex (P) falls below the breakup threshold, cf. 
similar instability scenarios reported in Refs. [8]- 12]. 

In region E separating the stability and breakup do- 
mains in Fig. 1, P exceeds the breakup threshold, hence 
the vortex does not split. However, in this region the 
evolution leads to squeezing the circular vortex into an 
elliptic rotating one, which remains stable. For instance, 
in Fig. 3, the circular vortex persists until z — 3700. The 
further evolution, lasting for Az sa 200, ends up with the 
transition into a stable elliptic soliton of a larger ampli- 
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FIG. 3: (Color online) The self-trapping of a stably rotating 
elliptic vortex. The configurations displayed at z = 29968 and 
z — 30000 correspond to the rotation of the ellipse by 180°. 
Images in the bottom row display the established "volcanic" 
shape of the elliptic vortex. 



tude, which revolves with a constant period, Z « 32. 
Its profile features a "volcanic" shape, with an undu- 
late "crater" (see the bottom row). This is a novel type 
of vorticity-carrying solitons in dissipative media. Note 
that elliptically deformed vortices were observed in a con- 
servative nonlocal optical medium |13|, a. nd they were 
studied theoretically in diverse settings [lj| , but in those 
cases the ellipticity was imposed by anisotropic bound- 
aries, while here the transition to the elliptic shape is 
generated intrinsically. Rotating solitons of other types 
("dipole propellers") were also predicted in conservative 
media, but in those cases the rotation was imparted, 
rather than spontaneous |15| . 

In the right-bottom ("unstable") region of Fig. 1, the 
vortices suffer destruction after the propagation distance 
measured in hundreds of units (closer to the border of 
region E, this distance extends to thousands). Four re- 
gions D,C,B, and S, which separate the stable and un- 
stable regions in Fig. 1, feature other remarkable vortical 
patterns. In region D, the stable elliptic vortex sponta- 
neously becomes eccentric. As shown in Fig. 4, it devel- 
ops the eccentricity after z — 8500, and after z — 14000 
it slips out from the central position and starts precessing 
around it, thus featuring double rotation (the precession 
and spinning) in the same direction, with an apparently 
locked period ratio, 216/54 = 4, which may be caused 
by a nonlinear parametric resonance in the original spin- 
ning elliptic vortex. The eccentric elliptic vortices remain 
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FIG. 4: (Color online) A stable eccentric spinning vortex, 
which orbits around the center. 
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FIG. 5: (Color online) The evolution of a stable rotating 
crescent-shaped vortex. 



stable in the course of this complex motion. 

Another variety of spontaneously deformed rotating 
vortices occupies region C in Fig. 1. As shown in Fig. 
5, after z — 660 the circular symmetry is broken by gen- 
erating a crescent-shaped vortex, which fills a half of the 
original circle. The crescent rotates around the center, 
with period Z k, \2 (see panels in Fig. 5 for z = 25988 
and z — 26000). Stable crescent-shaped solitons were 
predicted in the conservative model of rotating BECs 
trapped in anharmonic axisymmetric potentials |16| . 

The above species of the vortices feature stable shapes. 
On the other hand, in region B of Fig. 1 they are 
subject to an oscillatory instability, which immediately 
transforms them into robust breathers, which keep the 
vorticity and axial symmetry. Such stable breathers are 
similar to the one shown in the first panel of Fig. 6 (the 
period of its oscillations is Z ~ 5). Breathing vortices 
also appear in region S*; however, as shown in Fig. 6, 
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FIG. 6: (Color online) A vortex breather formed from the VA- 
predicted input at z ~ 20 and shown from above a,t z — 400, 
evolves into an elliptically-shaped "slanted volcano", whose 
profile is displayed at 2: = 25988. It rotates with period 
Z « 26 (the images shown at z = 25974 and z = 26000 
are separated by the period). 



they spontaneously develop the elliptic deformation and 
cease intrinsic oscillations at z ~ 500. The resulting vor- 
tex develops into an elliptic "volcano" which rotates with 
period Z « 26. The difference from Fig. 3 is that this 
vortex species is characterized by the "crater" which is 
not undulate but canted. 



Conclusion. Using the variational approximation and 
numerical simulations, we have developed the analysis 
of the general 2D model of optical media with the lo- 
calized linear gain compensated by losses. The balance 
of the diffraction, self-focusing, saturable nonlinear gain 
and the effective inhomogeneous loss provides for the sta- 
bilization of vortex solitons, which were unstable in previ- 
ously studied models of isotropic optical media. In addi- 
tion to the broad stability area for the vortices, stability 
domains are also found for several novel species of vortical 
solitons, including elliptic spinning ones, double-rotating 
eccentric vortices, revolving crescents, and breathers. 

The work suggests extensions in other directions. 
First, the stabilization may be also provided by the "ice- 
berg of gain" elevating above the surface of the "sea of 
loss", i.e., with the gain dominating near the center. In 
that case, it is sufficient to have the usual cubic loss, 
cf. the ID model developed in Ref. [6|]. On the other 
hand, for the stabilization of vortices it may be relevant 
to consider the loss-modulation profile with a minimum 
at a finite distance from the center. It is also interesting 
to consider a possibility of the stabilization of vortices by 
means of spatial modulations of the nonlinear loss, which 
may be relevant to optics and BEG alike [U^. Finally, 
a challenging problem is the stabilization of 3D vortex 
solitons. 
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